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Abstract 

The BABAR detector, which operates at the SLAC PEP-II asymmetric e^e^ cohider at energies 
near the T(45) resonance has cohected about 23M BB pairs in year 2000. Based on this data 
sample, we present a first study of sin2/3, with samples of B^ J/ipKg , B^ — > ip(2S)Kg and 
B° J/ipK^ decays. The measured value is sin2/3 = 0.34±0.20 (stat) ±0.05 (syst). In addition, we 
present preliminary measurements of charged and neutral B meson lifetimes and B^B^ oscillation 
frequency. 
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1 Introduction 



The primary goal of tiie BABAR experiment at PEP-II is to over-constrain the Unitarity Triangle. 
The sides of this triangle can be measured through non-CP violating physicSj_such as Vub, Vcb-, ^td 
measurements, while its angles are accessible through CP violating processeslil. 

2 PEP-II 

The PEP-II B Factory! is an e^e~ colliding beam storage ring complex on the SLAC site designed to 
produce a luminosity of at least 3x10^^ cm^^s^^ at a center ~of~mass energy of 10.58 GeV, the mass 
of the T(45') resonance. In the 2000 run, the achieved average luminosity was 3.3x10^^ cm~^s~^. 
The total collected luminosity was about 23fb~^. The machine is asymmetric with a High Energy 
Ring (HER) for the 9.0 GeV electron beam and a Low Energy Ring (LER) for the 3.1 GeV positron 
beam. This corresponds to /37=0.56 and makes it possible to measure time dependent CP violating 
asymmetries. It corresponds to an average separation of (3^cT=2bQ /xm between the two B mesons 
vertices. 

3 BABAR 

3.1 Detector description 

The BABAR detector is described ini. The volume within the BABAR superconducting solenoid, 
which produces a 1.5 T axial magnetic field, consists of: a five layer silicon strip vertex detector 
(SVT), a central drift chamber (DCH), a quartz-bar Cherenkov radiation detector (DIRC) and a 
Csl crystal electromagnetic calorimeter (EMC). Two layers of cylindrical resistive plate counters 
(RPCs) are located between the barrel calorimeter and the magnet cryostat. All the detectors 
located inside the magnet have full acceptance in azimuth. The integrated flux return (IFR) 
outside the cryostat is composed of 18 layers of steel, which successively increase in thickness away 
from the interaction point, and are instrumented with 19 layers of planar RPCs in the barrel and 
18 in the endcaps. 

3.2 Event reconstruction 

Charged particles are detected and their momentum is measured by a combination of the DCH 
and SVT. The charged particle momentum resolution is approximately given by {Spt/pt)"^ = 
(0.0015pr)^ + (0.005)^, where px is in GeV/c. The SVT, with a typical resolution of 10 /im per hit, 
provides excellent vertex resolution both in the transverse plane and in z. The vertex resolution 
in z is typically 50 fim for a fully reconstructed B meson and of order 100 for the distance 
among the two B mesons when only one is fully reconstructed. Leptons and hadrons are identified 
using a combination of measurements from all the BABAR components, including the energy loss 
dE/dx in the helium-based gas of the DCH (40 samples maximum) and in the silicon of the SVT 
(5 samples maximum). Electrons and photons are identified in the barrel and the forward regions 
by the EMC, and muons are identified in the IFR. In the barrel region the DIRC provides excellent 
kaon identification over the full momentum range above 250MeV/c. 



4 sin2/3 measurement 



In e"^e storage rings operating at the T(4S') resonance a pair produced in a T(4S') decay 

evolves in a coherent P-wave until one of the B mesons decays. If one of the B mesons {Btag) 
can be ascertained to decay to a state of known flavor at a certain time ttag, the other B (Bcp) 
is at that time known to be of the opposite flavor. For the measurement of sin2/3, Bcp is fully 
reconstructed in a CP eigenstate {J/t/jK^, i/j{2S)Kg or J/iJjK^). By measuring the proper time 
interval At = tcp — ttag from the Btag decay time to the decay of the Bqp (tcp), it is possible to 
determine the time evolution of the initially pure B^ or B^ state: 

f±{At; r, Amd, 2? sin 2/3) = ire"^!^*! [1 =F Vr]cpsm2[3 x smArudAt] , (1) 

where the + or — sign indicates whether the Btag is tagged as a, B^ oi & B^ , respectively. The 
dilution factor T> is given by I? = 1 — 2w, where w is the mistag fraction, i.e., the probability that 
the flavor of the tagging B is identified incorrectly, r/cp is the CP eigenstate of the final state and 
it is r/cp = —1 for the J/ipK^ and tlj{2S)Kg modes, rjcp = +1 for the J/ipK^ mode. Although 
less pure, the J/ijjK^ mode is very important because the oscillation is expected to be opposite 
to the other ones. A direct CP violation term proportional to cos Arud At could arise from the 
interference between two decay mechanisms with different weak phases. In the Standard Model, 
we consider that the dominant diagrams for the decay modes have no relative weak phase, so no 
such term is expected. 

To account for the finite resolution of the detector, the time-dependent distributions f± for B^ 
and B^ tagged events (Eq. ||) must be convoluted with a time resolution function TZ{At;a): 

J^±( At ; r, Anid, Vrjcp sin 2/3, a ) = /±( At ; T, Airid, Vrjcp sin 2/3 ) ® 7^( At ; a ) , (2) 

where a represents the set of parameters that describe the resolution function. 

Finally, the time-dependent distributions need to account for the background with additional 
parameters that characterize both its sample composition and the At distribution. The rjcp = — 1 
background components are parametrized with an ARGUS function and a small peaking component 
in the energy substituted mass mES- While the ARGUS component is extracted from data, the 
peaking one is derived from data. The rjcp = +1 background is made of a component from B decays 
in J/ip , whose parameters are taken from MC, and the other sources which are characterized using 
the J/tp mass sidebands. 

Since no time-integrated CP asymmetry effect is expected, an analysis of the time-dependent 
asymmetry is necessary. 

4-1 Analysis 

For this analysis, published in§, we use a sample of 23fb~^ of data recorded in year 2000, of which 
2.6 fb~^ was recorded 40MeV below the T{4S) resonance (off-resonance data). 
The measurement of the CP-violating asymmetry has five main components : 

• Selection of the signal ^/B° ^ J/4'K'i, B^/B^ i^{2S)K^ and P°/-B° ^ J/tp events, 
as described in detail inQ. 

The selection of the two modes with the Kg can profit from the fact that there are two dis- 
criminating kinematic variables, AE, the difference between the reconstructed and expected 




Figure 1: J/ipK^ (K^ tt+tt-, tt'^tt"; and tpK^ (K'^ tt°tt°) signal (top) and J/tpKl signal (bottom) . 



B meson energy measured in the center-of-mass frame, and mES) the beam-energy substi- 
tuted mass (see Fig. |l|). The background is mainly coming from combinatorics in continuum 
and other B decays and its properties can be estimated from the m-Es sidebands. 

The J/V' i^^sample is instead less pure because the the momentum is not reconstructed. 
The B mass constraint is therefore imposed and only one discriminating variable is left ( /S.E 
) as shown in Fig. |^. The background is mainly due to other B J/tpX decays and its 
properties are determined on MC. 

Signal event yields and purities, determined from a fit to the m-Es distributions after selection 
on ^E, are summarized in Table 

• Selection of decays in flavor eigenstate (i?flav)- candidates are formed by combining a D* 
or with a tt'^ , p'^ (p^ ''^^) ^ (a^ 7r'^7r~7r'''), or by combining a J/ip candidate 
with a K*^ {K*^ K~TT^p. Their background is mainly due to conbinatorics and can be 
studied in the ttt-es sidebands. Yields and purities are also summarized in Table ||. 

• Measurement of the distance Az between the vertex of the reconstructed B meson (i?rec) and 
the vertex of the flavor-tagging B meson (Stag)- 

In the reconstruction of the B^^c vertex, we use all charged daughter tracks. The vertex for 
the -Stag decay is constructed from all the remaining tracks in the event. 

In order to reduce bias and tails due to long-lived particles, and candidates are used 
as input to the fit in place of their daughters. In addition, tracks consistent with photon 
conversions (7 — > e~^e~) are excluded. To reduce contributions from charm decay products, 
which bias the determination of the vertex position, the track with the largest vertex con- 
tribution greater than 6 is removed and the fit is redone until no track fails the requirement 
or only one track remains. 



Table 1: Number of tagged events, signal purity and result of fitting for CP asymmetries in the full CP sample and 
in various subsamples, as well as in the Bflav and charged B control samples. 



Odliipit! 


A/1 


r yxiiiy \/Q) 






273 


96 ± 1 


0.25±0.22 


Jli\)Kl 


256 


39 ± 6 


0.87±0.51 


Full CP sample 


529 


69 ± 2 


0.34±0.20 


Jl^K'l, i){2S)Kl only 








Jli^Kg (Kg TT^ir ) 


188 


98 ± 1 


0.25±0.26 


J/ipKg [Kg TT IT ) 
ij{2S)K°g [K^g TT+TT-) 


A 1 

41 


or 1 

85 ± b 


-U.L)5±U.b6 


44 


97 ±3 


0.40±0.50 


Lepton tags 


34 


99 ± 2 


0.07±0.43 


Kaon tags 


156 


96 ± 2 


0.40±0.29 


NTl tags 


28 


97 ±3 


-0.03±0.67 


MT2 tags 


55 


96 ±3 


0.09±0.76 


-Bflav sample 


4637 


86 ± 1 


0.03±0.05 


Charged B sample 


5165 


90 ± 1 


0.02±0.05 



Prom the measurement of Az and of the B momentum, At can be computed. At an 
asymmetric-energy B Factory, in fact, the proper decay-time difference At is, to an excellent 
approximation, proportional to the distance Az between the two i?'^-decay vertices along the 
axis of the boost. At ~ Az/c (/?7). 

The time resolution function in equation ^ is described by a sum of three Gaussian distribu- 
tions (called the core, tail and outlier components) with different means and widths: 

For the core and tail Gaussians, the widths are scaled by the event-by-event measurement 
error ciAt derived from the vertex fits: a\^2 = Sx^^ x fAi- In data, approximately 65% of the 
area of the resolution function is in the core Gaussian. The width of the core Gaussian is 
approximately 110 /im or 0.7 ps; the width of the tail Gaussian is approximately 300 /im or 
1.8 ps. The third Gaussian has a fixed width of 8 ps and no offset; it accounts for the fewer 
than 1% of events with incorrectly reconstructed vertices. 

• Determination of the flavor of the Bt^g ■ 

Each event with a CP candidate is assigned a B^ or B'^ tag if the rest of the event (i.e., with 
the daughter tracks of the Bcp removed) satisfies the criteria from one of several tagging 
categories. In other words, a B'^ tag indicates that the Bqp candidate was in a B'^ state at 
At = 0; a B^ tag indicates that the Bcp candidate was in a S° state. 

Two tagging categories rely on the presence of a fast lepton (Lepton category) and/or one 
or more charged kaons in the event (Kaon category). Two categories, called neural network 
categories (NTl and NT2), are based upon the output value of a neural network algorithm 
applied to events that have not already been assigned to lepton or kaon tagging categories. 



Table 2: Mistag fractions measured from a maximum-likelihood fit to the time distribution for the fully-reconstructed 

B° sample. The uncertainties on e and Q are statistical only. 



Category 


e{%) 


w (%) 


Aw (%) 


Q{%) 


Lepton 


10.9 ±0.4 


11.6 ±2.0 


3.1 ±3.1 


6.4 ±0.7 


Kaon 


36.5 ±0.7 


17.1 ± 1.3 


-1.9± 1.9 


15.8 ± 1.3 


NTl 


7.7 ±0.4 


21.2 ±2.9 


7.8 ±4.2 


2.6 ±0.5 


NT2 


13.7 ±0.5 


31.7 ± 2.6 


-4.7 ±3.5 


1.8 ±0.5 


All 


68.9 ± 1.0 






26.7 ± 1.6 



The figure of merit for each tagging category is the effective tagging efficiency Qi = Si (1 — 2wif' , 
where is the fraction of events assigned to category i and Wi is the mistag fraction. The 



effective tagging efficiency as evaluated in data is summarized in Tab. 4.1 



The mistag fractions and the tagging efficiencies obtained by combining the results from 
maximum likelihood fits to the time distributions in the hadronic and semileptonic samples 



are summarized in Table 4.1 



Extraction of the amplitude of the CP asymmetry and the value of sin2/3 with an unbinned 
maximum likelihood fit. In order to extract as much information from the data itself and 
properly account for correlation, the fit is performed simultaneously to the CP and the flavor 
eigenstates. There are 35 parameters free in the fit: 

— Value of sin2/3; 

— Signal resolution function: Nine parameters Oj to describe the resolution function 
for the signal, being scale factors Si^2 for the event-by-event Az resolution errors of the 
core and tail Gaussian components, individual core biases Si^i per tagging category and 
a common tail bias 82, and the tail /i and outlier fractions; the width of the outlier 
component is taken to be a fixed 8 ps with zero bias; 

— Signal dilutions: Eight parameters to describe the measured average dilutions (Pj) 
and dilution differences 5T>i in each tagging category. 

— Background resolution function: Three parameters are used to describe a common 
resolution function for all non-peaking backgrounds, which is taken as a single Gaussian 
distribution with a scale factor Si for the event-by-event Az errors and an common bias 
(^1, and an outlier fraction /s; the width of the outlier component is taken to be a fixed 
8 ps with zero bias; 

— -Bfiav background properties: A total of 13 parameters describe the -Bfiav background 
properties. We make several assumptions to simplify the parameterization of the back- 
ground contributions and assign a corresponding systematic uncertainty. The mixing 
background contribution is assumed to be absent, = 0. The size of the peaking 
background is determined from Monte Carlo simulation to be (^pc^k ~ -'-•^ ^ 0.5% of the 
signal contribution in each tagging category. This contribution is dominantely from 
events, so Am^ = 0, E^^p^j^j^ = E^+ and -D^peak taken from the B^ data sample. The 
effective dilutions for the prompt (-Dff^, 4 parameters) and lifetime (i^ff^, 4 parame- 
ters) contributions are allowed to vary. The relative amount of these two contributions 
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Figure 2; Variation of the log likelihood as a function o/sin2/3 (left), for the whole sample (full curve), the rjcp = —1 
sample (dashed curve) and the rjcp = 1 sample (dotted line). Distribution of At for (a) the B" tagged events and (b) 

the B'^ tagged events in the CP sample (right). 



is allowed to vary, independently in each tagging category (4 parameters). For the life- 
time contribution, r?|^ is assumed to be same for all tagging categories, giving one free 
parameter. 

— CP background properties: One parameter, the fraction of prompt relative to life- 
time background, assumed to be the same for each tagging category, is allowed to float 
to describe the CP background properties. The effective dilutions of the lifetime and 
peaking contribution are set to zero (-Dj^f = Df^peak ~ O), corresponding to no CP- 
asymmetry in the background. The size and parameters of the peaking background is 
again determined from Monte Carlo simulation. The fraction of peaking background is 
'^peak = 1 i 1% of the signal contribution, independent of tagging category. This contri- 
bution is assumed to have dilutions and lifetime parameters in common with the signal 
contribution. Finally, the lifetime of the lifetime background is assumed to be r^o in all 
tagging categories. 

The maximum-likelihood fit for sin2/3, using the full tagged sample, gives: 

sin2/3 = 0.34 ± 0.20 (stat) ± 0.05 (syst) . (4) 

For this result, the lifetime and Ama are fixed to the current best values i. The log 
likelihood is shown as a function of sin2/3 and the At distributions for B'^ and B^ tags are 
shown in Fig. |2|. 

The dominant sources of systematic error are the assumed parameterization of the At resolution 
function (0.04), due in part to residual uncertainties in the SVT alignment, and uncertainties in 
the level, composition, and CP asymmetry of the background in the selected CP events (0.02). 
The systematic errors from uncertainties in Am^o and r^o and from the parameterization of the 



background in the selected Sflav sample are found to be negligible. An increase of 0.02 /ips~^ in 
the assumed value for Am^o decreases sin2/3 by 0.012. 

The large sample of reconstructed events allows a number of consistency checks, including 
separation of the data by decay mode, tagging category and i?tag flavor. The results of fits to these 
subsamples are shown in Table || for the high-purity events. Table |l] also shows results of fits 
with the samples of non-CP decay modes, where no statistically significant asymmetry is found. 

5 Measurements of charged and neutral B meson lifetimes and B^B^ oscillations 

These measurements can be used to test theoretical models of heavy-quark decays and to constrain 
the Unitarity Triangle (via the sensitivity to the value of the CKM matrix element Vtd)- 
One B {B 

rec) is fully reconstructed in an all-hadronic (5° D^*^ 7r+, D^*) p+, D^*^ af, 
J/ipK*'^ and B~^ — > D^*^^tt~^, J/\IjK~^, ijj(2S)K~^) modes. The number of selected events and 
purities are summarized in Table |l[ 

The measurement of At and, when needed, the tagging of the recoiling B is done in the same 
way as for the sin2/? measurement, so that these measurements constitute also a valuable validation 
of the CP measurement. 

5.1 Lifetime Measurements 

The B^ and B~^ lifetimes are extracted from a simultaneous unbinned maximum likelihood fit to 
the At distributions of the signal candidates, assuming a common resolution function. An empir- 
ical description of the At background shape is assumed, using niES sidebands with independent 
parameters for neutral and charged mesons. Fig. |3| shows the At distributions with the fit result 
superimposed. 
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Figure 3: At distributions for B'^ /B'^ (Isft) <^r,d B'^ /B (right) candidates in the signal region (mES> 5.27GeV/c^;. 
The result of the lifetime fit is superimposed. The background is shown by the hatched area. 



5.2 Time-dependent B^B^ mixing 

A time-dependent B^B^ mixing measurement requires the determination of the flavor of both B 
mesons. Considering the B^B^ system as a whole, one can classify the tagged events as mixed or 
unmixed depending on whether the Btag is tagged with the same flavor as the B^ec or with the 
opposite flavor. 



From the time-dependent rate of mixed (Nmix) and unmixed {Nunmix) events, the mixing 
asymmetry a(At) = [Nunmix — N mix) / unmix + ^mix) is calculated as a function of At and fit to 
the expected cosine distribution. A likelihood fit with 34 free parameters is performed. The free 
parameters are the same as in the sin2/3 fit, apart from sin2/? itself and the background to the CP 
eigenstates which are not used, and Am^, which is, of course, floated. 
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Figure 4: At distribution for mixed and unmixed events (left) and time- dependent asymmetry a[A.t) between unmixed 

and mixed events (rigth) 



Fig. I shows the At and a(At) = {Nunmix - Nmix)/{Nu 
result superimposed. 
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5. 3 Results 

The preliminary results for the B meson lifetimes are r^o = 1.546ib0.032 (stat) ±0.022 (syst) ps and 
tb+ = 1.673 ±0.032 (stat) ±0.022 (syst) ps and for their ratio is tb+/tbo = 1.082 ±0.026 (stat) ± 
0.011 (syst). 

We measure the B^B^ oscillation frequency: Am^ = 0.519 ± 0.020 (stat) ± 0.016(syst) ^ps"^ 
The above results are consistent with previous measurementsB and are of similar precision. The 
mixing result is compatible with a BABAR measurement using di-leptonsu. 

6 Conclusions 

We have presented BABAR's first measurement of the CP-violating asymmetry parameter sin2/3 in 
the B meson system: 

sin2/3 = 0.34 ± 0.20 (stat) ± 0.05 (syst) . (5) 

Our measurement is consistent with the world average sin2/3 = 0.9 ± 0.40, and is currently limited 
by the size of the CP sample. 

We have also presented time-dependent mixing and lifetime measurements, performed for the 
first time at the T(4S'). 
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